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We consider bosonic atoms in an optical lattice at integer filling, tuned to the superfiuid-Mott 
insulator critical point, and coupled to a single, mobile impurity atom of a different species. This 
setup is inspired by current experiments with quantum gas microscopes, which enable tracking 
of the impurity motion. We describe the evolution of the impurity motion from quantum wave 
packet spread at short times, to Brownian diffusion at long times. This dynamics is controlled by 
the interplay between dangerously irrelevant perturbations at the strongly-interacting field theory 
describing the superfluid-insulator transition in two spatial dimensions. 



I. INTRODUCTION 

The problem of a single particle interacting with its en- 
vironment is encountered frequently in condensed matter 
physics, and in many variations. Such scenarios appear 
in the form of genuine defects in materials, like mag- 
netic impurities in metals and the associated ubiquitous 
Kondo problem.^ Impurity problems are also often used 
as a stepping stone to gain insight into the behavior of 
strongly correlated systems, such as studying a single 
hole doped into a Mott insulator helps understanding 
some properties of high-temperature superconductors. 

In recent years, impurity problems have also started 
to gain some attention in the context of ultracold atomic 
gases. Amongst others, the fermionic polaron problem 
has been studied thoroughly in connection with strongly 
imbalanced Fermi gases close to a Feshbach resonance^"^, 
where a single "spin-down" atom interacting with a Fermi 
sea of "spin-up" atoms appears as a polaronic quasi- 
particle and eventually forms a two-particle bound-state 
with one of the majority atoms if the interactions are 
strong enough^"^^. More generally, such polaronic prob- 
lems appear when we consider the nature of the threshold 
of spectral functions of gapped excitations above many- 
body ground states. ^^'^"^ 

Here we are interested in a seemingly similar prob- 
lem, where a mobile impurity is coupled to a continuum 
of gapless excitations at a quantum critical point, in- 
stead of gapless particle-hole excitations at a Fermi sur- 
face. In particular, we want to study the situation where 
bosons in an optical lattice are tuned to the superfluid- 
Mott insulator transition, and a single atom of a different 
species is coupled to the bosons via a density-density in- 
teraction. This scenario is especially interesting due to 
the recent experimental implementation of quantum-gas 
microscopes^*'^^, which allow tracking the motion of the 
impurity through the bosonic bath. Indeed, a recent com- 
bined experimental and theoretical work has addressed a 
very similar situation in an effectively one-dimensional 
system^^ . We are going to focus on the other experimen- 
tally relevant, two dimensional situation instead. From 
our perspective, this case also happens to be the most 
interesting scenario, because in d = 2 spatial dimensions 



the critical behavior at the SF-MI transition is described 
by a non-trivial fixed point without well-defined quasi- 
particle excitations. 

The main motivation for our work is the prospect of 
studying aspects of quantum critical transport in a clean, 
well-defined experimental setting. Calculating transport 
coefficients of systems close to a quantum phase transi- 
tion remains a challenging problem in theoretical con- 
densed matter physics. The main reason is that all 
traditional methods for calculating transport properties 
rely on a quasiparticle picture of the low-energy excita- 
tions above the ground-state, which breaks down at a 
quantum critical point. Indeed, strongly coupled confor- 
mal field theories (CFT) that describe physical proper- 
ties of physical systems close to quantum critical points 
typically have no well defined quasiparticle excitations. 
Some progress has been made recently using the so called 
AdS/CFT correspondence, which circumvents this prob- 
lem by mapping the strongly coupled CFT to a dual, 
weakly coupled gravitational theory, where transport co- 
efficients can be computed reliably. 

In this paper our main interest is to calculate the dif- 
fusion constant D of the impurity as a function of tem- 
perature in the quantum critical regime, shown in Fig. 1. 
This diffusion constant can be measured directly in ex- 
periments with quantum-gas microscopes by evaluating 
the mean square displacement (x^) of the impurity atom 
from its initial position after a time t 

{^^)^4Dt, (1.1) 

where the brackets denote an average over many experi- 
mental realizations and the time t ^ t has to be larger 
than the typical collision time r, below which the im- 
purity propagates ballistically. Several theoretical works 
have addressed similar questions using the holographic 
correspondence in the context of a heavy quark mov- 
ing through a quark gluon plasma. ^^"^"^ Here we will use 
a conventional method based on the Boltzmann equa- 
tion, however, which turns out to provide an adequate 
description at low temperatures. The quantum critical- 
ity of the bulk critical point has a single characteristic 
time which controls transport and relaxational processes 
close to equilibrium tu ~ h/{kBT).'^^ We will show here 
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FIG. 1: Schematic phase diagram of the superfluid - Mott 
insulator transition of a Bose-Hubbard model in d = 2 spatial 
dimensions at integer filhng, as a function of temperature T 
and the tuning parameter r = — (see Eq. (2.10)). Tkt 
marks the Kosterhtz-Thouless transition temperature to the 
superfluid state. The quantum critical point, marked by XY, 
is in the universality class of the 3d XY model. We focus 
on transport properties of an impurity atom coupled to the 
bosons in the quantum critical regime. 



that the situation with impurity motion is not as uni- 
versal, and the impurity time scale is determined by an 
interplay between several couplings which are formally 
irrelevant at the critical point. One of the these cou- 
plings is the inverse effective mass, 1/m, of the mobile 
impurity; there is the associated energy scale mc^ , where 
c the velocity of the bulk bosonic excitations at the crit- 
ical point. The value of m controls the quantum wave 
packet spread of the impurity at short times. The other 
"irrelevant" parameters couple the impurity to the bulk 
excitations, and associated energy scales determine the 
diffusion of the impurity at long times even at the lowest 
temperatures. One of our interesting results is that the 
diffusion constant of the impurity does not depend on m 
at low temperatures, however. 

The rest of the paper is organized as follows. In sec- 
tion II we derive an effective low energy theory for a mo- 
bile impurity coupled to bosons at the superfluid-Mott- 
insulator (SF-MI) transition and we show that the pres- 
ence of three-body interactions changes the theory signifi- 
cantly at the particle-hole symmetric point. In section III 
we study perturbative self-energy corrections to the im- 
purity propagator and calculate its diffusion constant as 
a function of temperature, neglecting three-body interac- 
tions. We repeat these calculations including three-body 
interactions in section IV and show that this changes 
the results qualitatively at low temperatures. Finally we 
perform a renormalization group analysis of the effective 
field theory in Sec. V. 



II. FIELD THEORY APPROACH 

We start by deriving the effective low-energy theory 
for a Bose-Hubbard model at the SF-MI critical point 



in d = 2 spatial dimensions, following Ref. 26. A con- 
venient starting point is to consider particle- and hole- 
excitations on top of a Mott insulating state, which are 
described by the Lagrangian 
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(2.1) 



Here the fields p = p(x, r) and h = /i(x, t) represent 
particle- and hole excitations on top of a state with an 
integer filling of bosons. The respective energy gaps for 
these excitations are denoted by Ap and A^, their effec- 
tive masses by rup and ruh- The last term in Eq. (2.1) 
creates or annihilates particles and holes in pairs, as re- 
quired by the conservation of the total number of bosons, 
and we do not explicitly show higher order terms in the 
fields p and h that are allowed by symmetry. 

Ultimately we are interested in the experimentally rele- 
vant situation where the transition happens at a constant 
density, i.e. at the tip of the Mott-lobe. This happens 
only if the gaps Ap and A/j are equal, such that parti- 
cles and holes condense at the same time when the gap 
closes. For the moment we will consider the general case, 
however, and parametrize the gaps as 



Ap = A -5/2 
A,, = A + 5/2 



(2.2) 
(2.3) 



Now we add a single, mobile impurity atom to the bosons, 
the dynamics of which is described by the free particle 
Lagrangian 



(2.4) 



where the field 6(x, r) represents the impurity and the 
chemical potential < has to be adjusted such that 
the density of impurities is zero. This is important only 
at finite temperatures, however, where we need to ensure 
that no artificial, thermally excited impurities exist. At 
zero temperature we can safely set /i = 0. 

Within a microscopic Bose-Hubbard like model the 
coupling between the impurity and the bosons is an onsite 
density-density interaction, which is a pure two-body in- 
teraction in the simplest case. Since the local boson den- 
sity is given by the difference between particle- and hole- 
densities (up to the constant average number of bosons 
that we neglect in the following, as it only renormalizes 
the impurity's chemical potential), the interaction term 
in our effective field theory takes the form 



,(2) 
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(2.5) 



where the subscript 2 of the coupling constant U2 in- 
dicates that this term derives from a two-body interac- 
tion. Furthermore we include a term that descends from 
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a three-body interaction 

c\^=u,\b\'{\p\'^\h\Y , (2.6) 

because in the case of particle-hole symmetry (S = 0) this 
interaction generates couplings that are more relevant in 
an RG sense than those deriving from two-body interac- 
tions only, as will be shown below. To avoid confusion we 
note that a three-body interaction in the Bose-Hubbard 
model of the form ^ b\hi ni(ni — 1), where Ui denotes the 
number of bosons on site i and b\ is the impurity cre- 
ation operator, also generates a term of the form (2.5) in 
addition to (2.6). For brevity we will always refer to the 
term in (2.6) as the three-body interaction, however. 

Eventually our system is described by the total La- 
grangian 

C — Cc + -^imp + ■Cjjjt -|- . (2-7) 

In order to arrive at the more familiar low energy theory 
for the SF-MI critical point we define the fields 

i, = {p + h*)/V2 (2.8) 
C - {p-h*)/V2 (2.9) 

and integrate out the field ^. After a rescaling of tp and 
defining r = — as well as = {A + X){mp + 
rufi) / (Ampmii) we arrive at our final result 

£ = ^P* (S dr - - c^V^ + r^iP + ^\iP\* 

+b*(dr-^-^i)b+iSu2+U3)\b\''\4,\^ 
\ 2m / 

+U2{'4>drij* -^p*dr^)\bf . (2.10) 

The terms shown are the leading order terms within a 
gradient expansion. For the case of particle-hole sym- 
metry {S = 0, i.e. Ap — Afi) the dynamical critical ex- 
ponent of the bosons is z = 1 and the critical point is 
described by the well known C/(l) CFT in d -I- 1 dimen- 
sions. Even more interesting, the direct interaction term 
between the bosons and the impurity ~ Ibpl-^p vanishes 
if only two-body interactions are present. This somewhat 
counterintuitive result (in the sense that the particle-hole 
symmetric interaction term vanishes if the i/j^-theory is 
particle-hole symmetric) can be understood from the fact 

(2) 

that the two-body interaction term Cl^^l. is invariant un- 
der a particle-hole transformation {ip ^ ip* , (5 — > —S) 
only in combination with U2 — >■ —U2. This is not true for 
the three-body interaction term Cl^l, however, and that 
is why the parametrically smaller generates a term 
that is more relevant than the one generated by two-body 
interactions at the particle- hole symmetric point. The 
leading order terms that are consistent with the symme- 
try 

C[ip;S,U2,U3] = C[tp*;-S,-U2,U3] (2-11) 
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FIG. 2: Self energy contributions to the impurity propagator 
up to second order in the interaction U2. Solid lines denote 
the boson propagator, dashed lines correspond to the impurity 
propagator. 



are the ones shown in Eq. (2.10). Also note that the two- 
body interaction U2 couples the impurity density |6p to 
the time-like component of the bosonic current-density 

We stress again that we focus solely on the particle- 
hole symmetric case {S = 0) in the rest of the paper. A 
tree-level scaling analysis of the Lagrangian (2.10) using 
the exact critical exponents of the bosonic CFT shows 
that both interactions U2 and U3 are irrelevant in d = 2 
spatial dimensions at the p/h-symmetric point. Indeed, 
performing a scaling with z = 1 gives the scaling dimen- 
sions [ip] — (d — 1 + rj)/2 and \b] — d/2 for the fields, as 
well as [IV'P] — d + 1 ~ 1/v and [j] = d for the compos- 
ite bosonic density- and current operators per definition. 
Here 77 and v denote the anomalous dimension and the 
correlation length exponent of the bosonic CFT, which 
take the values 77 w 0.038 and v 0.672 at the quantum 
critical point in d = 2 spatial dimensions. From these 
scaling dimensions we get [U2] = l — d and [u^] = l/v — d, 
thus both couplings are irrelevant in d = 2. Also note 
that the curvature of the impurity dispersion is an irrele- 
vant perturbation, as the scaling dimension of the impu- 
rity mass m is [m] = 1. We will perform a more elaborate 
RG analysis in Sec. V. 



III. TWO-BODY INTERACTIONS 

In this section we analyze the field theory (2.10) at 
the particle-hole symmetric point [S — 0), assuming that 
three body interactions are absent (i.e. U3 — 0). The 
Lagrangian (2.10) thus takes the simplified form 

£ = r{-d^-c^^^+r'^^+^\^Pf (3.1) 

+b* {9r-^)b + U2{^ - rdr^') l&l' • 

Note that we dropped the chemical potential since it is 
not important at zero temperature, as argued above. 

A. Perturbative analysis (T=0) 

As a first step towards an understanding of the im- 
purity's properties we calculate its self-energy perturba- 
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tively in the interaction strength U2 at zero temperature 
in the following. Since U2 is irrelevant, we expect the per- 
turbative calculation to give qualitatively correct results. 
The terms up to second order are shown in Fig. 2. The 
first order tadpole contribution turns out to be identi- 
cally zero. Note that the boson loop in the second order 
diagram corresponds to the time-like component of the 
bosonic current-current correlation function of the free 
bosonic theory due to the particular form of the interac- 
tion vertex in (3.1). At T = the form of the bosonic 
current-current correlator of the interacting i/j'^-theory at 
the critical point in d = 2 is fixed by conformal invariance 
and is given by 



^00 



(w + iO+)2 



(3.2) 



Here k is a numerical factor that is renormalized by the 
boson interaction g and takes the value k = 1/16 for the 



free bosonic theory. The functional form of Koo{q,uj) is 
independent of g, however. For this reason our results 
for the impurity self-energy below are perturbative in U2, 
but valid to arbitrary order in the boson interaction g. 

At zero temperature the imaginary part of the retarded 
self-energy up to second order in U2 is given by 

ImE_R(k, w) = / 0(w - £k-q)Imiir,fo(q>'^ ^ £k-q) , 

(3.3) 

where Ek = k'^ /{2m) denotes the impurity dispersion. At 
k = the integral can be evaluated exactly and we obtain 



ImE_R(0,a;) 



w + mc 



log 1 



2a; 



- LO 



37rc4 



(3.4) 



For finite momenta we can derive approximate analytic expressions, which take the form 

2uIk m4(29fc2 -I- 31m2c2^ 



ImSfl(k,M 



45^2 7r2 (m2c2-yfc2)3 
SujK (fc — mcY 



5uj^ + 0{5uj'^) for k <mc 



I 3^2 7r2 



with (5w = w — Ek > 0. 



(3.5) 



+ 0{SUJ, (fc - TOC)") 



for k > mc 



r 



The qualitative change of the self-energy at A: = mc can 
be understood from kinematic considerations. The only 
lifetime limiting process for the impurity at zero temper- 
ature happens to be the excitation of a bosonic " phonon" 
mode. Due to energy- and momentum conservation this 
process is prohibited as long as the group velocity of the 
impurity Vg = 9fcek — k/m is smaller than the "sound" 
velocity c of the bosons. For this reason the impurity 
spectral function has a sharp delta-function peak for all 
momenta k < mc, which is reflected by the fact that the 
imaginary part of the self energy scales as ~ i5a;^ close to 
the quasiparticle pole. On the other hand, for k > mc 
the impurity can scatter strongly by exciting the bosons. 
Consequently the imaginary part takes a non-zero on- 
shell value, implying that there is no well defined im- 
purity excitation anymore. The appearance of this kine- 
matic constraint can be readily seen from Eq. (3.3), where 
the imaginary part of the current correlator is non-zero 
only if a; > £k-q + CQ- 

It is important to note that this threshold behavior for 
creating bosonic excitations at A: > mc holds to all orders 
in perturbation theory in the impurity-boson coupling 
U2- Higher order diagrams with more boson propagators 
describe processes where the impurity excites multiple 
boson modes, which in turn have an even higher ener- 
getic threshold. From a theoretical perspective it is an 



interesting question, if the impurity's properties above 
this kinematic threshold at k = mc are described by a 
RG fixed point. This is not the case, however, and we 
will come back to this problem in Sec. V. 



B. Diffusion constant from Boltzmann equation 

An important transport coefficient that can be directly 
measured in experiments is the diffusion constant D of 
the impurity. Here we calculate D{T,m, A) as a func- 
tion of temperature T, impurity mass m and the Mott 
gap A using the Boltzmann equation, which gives qual- 
itatively correct results as long as the impurity has a 
well defined quasiparticle peak. This condition is sat- 
isfied in the temperature regime ksT ^ mc^, where 
the kinematic threshold discussed in the previous sec- 
tion doesn't matter. Let us quickly estimate the experi- 
mentally relevant temperature regime before proceeding 
with the Boltzmann calculation. Typical temperatures 
in experiments^'*'^''' in units of the Hubbard interaction 
U are on the order of T/U « 0.1. Using the theoretical 
results for the hopping amplitude J as well as the sound 
velocity c at the critical point of the 2d Bosc-Hubbard 
model2^, which take the values (J/J7)crit — 0.0597 and 
(c/J)crit — 4.8, together with the band mass m — 2/J 
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of the impurity, we obtain / {m(?) « 0.15, which is 
typically well below the kinematic threshold. 

We are interested in a situation without external forces 
but with an initial density gradient, thus the Boltzmann 
equation takes the form 

5Jk(x, i) + - ■ Vx /k(x, t) = /[/k(x, t)] , (3.6) 
m 

with /k(x, i) as the Wigner distribution function of the 
impurity at momentum k, position x and time t. The 
collision integral has the standard form 



/[/k] = - ^ (l^k,k'/k - l^k'.k /k') , 



(3.7) 




2.0 0.0 



where the transition rates Wk,k' can be calculated using 
Fermi's golden rule and are given by 



Wk,i 



2 ul [lmifo« (k' - k, Sk' - ek) x 

xnB(ek' - Ek) 6(ek' - ^k) 
+Imii:oo(k' - k,£k - £k') X 

X (l + ni3(ek - Ek')) 0(£k - £k') • (3.' 



Here K^^{q^,ijj) denotes the time-like component of the 
retarded bosonic current-current correlation function, ub 
is the Bosc-Einstein distribution function and denotes 
the unit step function. Details of the derivation are given 
in appendix A. 

In typical experimental situations the impurity is ini- 
tially localized, released at some time t = Q and then is 
allowed to propagate through the bosonic bath. In prin- 
ciple the Boltzmann equation (3.6) allows to calculate 
the full time evolution of the impurity distribution func- 
tion. At short times collisions don't play a role and the 
impurity propagates ballistically 



/k(x, t < t) /k(x - kt/m, 0) 



(3.9) 



where r denotes the typical time between collisions. If we 
start aX t — with the particles localized in a Gaussian 
wavefunction 



FIG. 3: (Color online) Scaling function '^{^^, for the 

diffusion constant in Eq. (3.15). The scaling function has 
been evaluated numerically from Eq. (A12) using the one- 
loop expression for the current-current correlation function 
Koo{(\,u)) at finite temperature. 



takes the form (we set fc^ = 1 from now on) 

with a Maxwell-Boltzmann distribution of the momenta. 
Indeed, in the long time limit, where the Wigner distri- 
bution function of the impurity is close to equilibrium, 
we can derive a diffusion equation from the Boltzmann 
equation (3.6) for arbitrary linear collision integrals of the 
form (3.7). Again, details can be found in appendix A. 
The corresponding impurity diffusion constant D takes 
the simple form 



D = r 

m 



(3.13) 



with a transport relaxation time t which is given by 



k,k' 



k' 

1 - cos eik^k' 



/•Boltz 
/k 



(3.14) 



V'(x,0) 



1 



exp 



X 



(3.10) 



then the short time evolution of the Wigner distribution 
function according to Eq. (3.9) is free particle behavior 



/k(x,t«T) ^ j d? 



yg-^k-y^* (x--,t)i/;(x+-,t 



= 4cxp X 1] - k^a^ (3.11 



m 



At long times t ^ t collisions dominate and the impu- 
rity propagates diffusively. The Wigner distribution then 



where fc = |k| and ^ exp(— /3ek) is the equilibrium 

Boltzmann distribution function. Using the explicit form 
of the transition rates Wk,k' in Eq- (3-8) we can cast 
the expression for the diffusion constant D{T, m, A) in a 
scaling form 



D = 



T A 



(3.15) 



The corresponding scaling function $ can be calculated 
numerically and is shown in Fig. 3. Here we used the 
one loop result for the bosonic current correlator Kqq at 
finite temperature, where the dispersion relation of the 
bosonic modes = [ck)"^ -f acquires a tempera- 
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ture dependent mass term 



mt = 2riog 



V(2T) ^ V4 + e^/^ 



(3.16) 



which follows from a 1 /A'^ expansion in the limit — > oo 
of the bosonic CFT.^f 

At a fixed temperature T, the diffusion constant D 
takes its minimal value in the quantum critical regime 
directly above the quantum critical point, where the Mott 
gap vanishes (A = 0). This is consistent with the naive 
expectation that the phase space for scattering processes 
between the impurity and the bosonic modes is large, 
if the gap of the bosonic modes is small. The scaling 
function in the quantum critical regime for A = is 
shown in more detail in Fig. 4. Note that $ approaches 
a finite value <i>(0, 0) ~ 0.267 at zero temperature and 
scales linearly for small T/(mc^) ^ 1, as 



T 



,0 



0.267 + 4.4- 



T 



(3.17) 



Interestingly, the fact that the scaling function takes a 
constant value at T = and A = implies that the dif- 
fusion constant at the critical point does not depend on 
the impurity mass m at very low temperatures. Indeed, 
the curvature of the impurity dispersion, i.e. the inverse 
mass 7Ti~^, is an irrelevant perturbation, and here there is 
no conspiracy between the two irrelevant couplings 
and U2 that potentially could change the T dependence 
of D from that expected from its 1/^2 dependence. Only 
at higher temperatures the mass starts to play a role, 
leading to a crossover from a. D ^ behavior at very 
low temperatures to a. D ^ T^^ dependence at tempera- 
tures T > 0.06 mc^. Note that this crossover temperature 
is below the typical temperatures in experiments, where 
T > 0.15 mc^, thus we expect that the D ^ T^^ behavior 
should be observed. 

On a formal level, the fact that the scaling function at 
the critical point $(0,0) is constant can be traced back 
to the behavior of the current correlator at small frequen- 
cies, which scales as lmKQQ{q,uj) ~ w for uj <ti T. It is 
important to stress that this behavior, which is related 
to the fact that Imi4r^(q, a;) = — Imi4r^(q, —oj) is an odd 
function of frequency and expected to be an analytic at 
w = 0, is likely valid in general and not related to our 
one-loop approximation of the current correlator. Indeed 
it has been argued that the current-current correlation 
function takes a hydrodynamic form for uj,q ^ T 



<o(q,w) 



X 



Dq^ 



(3.18) 



the imaginary part of which obviously scales as ~ a; for 
small frequencies. This has also been confirmed by cal- 
culations using the AdS/CFT correspondence.^^ 

In the Mott regime away from the critical point, where 
the excitation gap A is finite, the scaling function di- 




T/mc^ 



FIG. 4: (Color online) Scaling function <E>(^^, O) for the dif- 
fusion constant in Eq. (3.15) as a function of temperature in 
the quantum critical regime above the critical point, where 
the Mott gap vanishes (A = 0). The scaling function has 
been evaluated numerically from Eq. (A12) using the one- 
loop expression for the current-current correlation function 
ifoo(q, i^) at finite temperature. The inset is an expanded 
view of the low temperature region, showing that the scaling 
function ${x,0) is non-zero at a; = 0. The dashed red and 
dash-dotted green lines show the scaling function in the Mott 
regime for A/{mc'^) — 0.05 and 0.1 as a function of T/{mc?), 
respectively. 



verges exponentially as the temperature approaches zero 



V 9 ' ') 



(3.19) 



This behavior can be understood from the fact that the 
impurity does not scatter with bosonic excitations at 
temperatures well below the Mott gap and thus prop- 
agates ballistically, which translates to a diverging diffu- 
sion constant. 

We also note that the mobility /i of the impurity is in- 
dependent of its mass at the critical point at low enough 
temperatures. The mobility is defined as /i = Vd/F, 
where Vd is the impurity's terminal drift velocity in re- 
sponse to an external force F. It is related to the diffusion 
constant D via the Einstein relation D = and thus 
has the same dependence on the impurity mass as the 
diffusion constant. 



IV. THREE-BODY INTERACTIONS 

We discussed in Sec. II that three-body interactions are 
more relevant than two-body interactions at the particle- 
hole symmetric quantum critical point. Indeed, the im- 
purity's behavior will be dominated by three-body in- 
teractions at low enough temperatures, albeit being 
parametrically smaller than U2. In this section we thus 
repeat the analysis of Sec. Ill for the case of three-body 
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interactions U3. For the sake of simplicity we set the two- 
body interaction to zero in this section and thus start our 
discussion from the Lagrangian in Eq. (2.10) with (5 = 
and U2 = 0. In reahty both interactions will be present, 
however, which will lead to a crossover in the tempera- 
ture dependence of the diffusion constant. The associ- 
ated crossover scale depends on the relative strength of 
the two interactions, which will be discussed briefly at 
the end of this section. 



B. Diffusion constant from Boltzmann equation 

For the computation of the impurity's diffusion con- 
stant the same arguments hold as in Sec. IIIB and we 
restrict our discussion here to the regime where the Mott 
gap vanishes (A = 0). The transition rates Wk,k' have 
the same form as in Eq. (3.8), with U2luiK^{q^,uj) re- 
placed by M§ Imnfl (q, w). Using the scaling form of the 
bosonic polarization function at finite temperature 



n(q, cj, T) = T^-^/" flicq/T, lu/T) 



(4.3) 



A. Perturbative analysis 



Again, we start by performing a perturbative analysis 
of the impurity's self-energy at zero temperature, includ- 
ing terms up to second order in . The two correspond- 
ing self-energy diagrams are analogous to the ones shown 
in Fig. 2. Note also that mixed terms ~ vanish by 
symmetry. We neglect the first order tadpole contribu- 
tion as it only renormalizes the chemical potential, which 
we set to zero in order to ensure a vanishing density of im- 
purities. The second order contribution to the self-energy 
can be computed by replacing the current-current corre- 
lator KoQ^q^uj) in Eq. (3.3) with the bosonic polarization 
hmction n(q, oj). Again we want to obtain a result for 
the impurity self-energy which is non-peturbative in the 
boson coupling g, thus we use the scaling form 



nfl(q,a;) 



122 
c q 



- ijJ' 



r 



(4.1) 



for the bosonic polarization function at zero tempera- 
ture, which follows from the scaling dimensions of the 
fields discussed in Sec. II. The corresponding integral 
to Eq. (3.3) for the imaginary part of the impurity self- 
energy can be evaluated exactly at k = 0. In the limit of 
small frequencies w <C mc^ we obtain 



ImEij(0, w) 



2 |sin(7r(3/2-lH)| 



47rc2(5/2- 1/i/) 



(4.2) 



Using the exact value for the correlation length 
exponent^^ v = 0.67155, the imaginary part of the self- 
energy scales as ~ o;^ "^^ at low frequencies, i.e. the im- 
purity spectral function again has a sharp delta-function 
peak. Note that this is not true in mean-field theory, 
where i^mf = 1/2 and the the imaginary part of the self- 
energy scales as ~ cj at small frequencies, which would re- 
sult in an algebraic singularity instead of a delta-function 
peak in the spectral function. 

At finite momenta we expect a similar behavior of ImS 
as described in Sec. Ill A for the case of two-body inter- 
actions, where the self-energy acquires a finite imaginary 
part beyond the kinematic threshold at A: = mc. We 
calculate this imaginary part explicitly in Sec. V. 



we can express the impurity diffusion constant in the scal- 
ing form 



D 



T 



(4.4) 



Unfortunately we cannot reliably calculate the scaling 
function $ in the equation above, because neither the 
mean-field approximation with v = 1/2, nor the leading 
order 1 jN approximation where v = \ are good approx- 
imations for the polarization function at finite temper- 
ature. In fact, from Eq. 4.1 one can see that nj^(q, w) 
is close to being constant at zero temperature, because 
the correlation length exponent almost takes the value 
V K, Iji at the critical point. However, it is reasonable 
to assume that the scaling function $(0;) in Eq. (4.4) 
will again take a constant value at a; = 0, because 
Imnij(q, cj) = — Imnij(q, — w) is an odd function of fre- 
quency and expected to scale linearly Imnj^(q, w) ^ to 
for small frequencies, similar to the case discussed in 
Sec. IIIB. Based on this assumption the diffusion con- 
stant scales as D ~ 7^-2.022 j^.^ temperatures. 

If both, two- and three-body interactions are present, 
there will be a crossover in the temperature dependence 
of the diffusion constant. Summing up the two inverse 
relaxation times using Matthiessen's rule, we obtain 



D- 



u\ T 



5-2/, 



"( 



1 + const. 



f2 rp2/v-2 

4 



(4.5) 



with a constant of order one. We thus expect a crossover 
from a T-2.022 rp-3 behavior at temperatures 

^ ^ (^i/"!)'^^'^"^'^'! provided that the mass-dependent 
crossovers happen at a higher temperature. 



V. RG ANALYSIS FOR THE THREE 
BODY-INTERACTION 

The simple tree level-scaling at the end of Sec. II 
showed that both U2 and U3 are irrelevant couplings at 
the particle-hole symmetric point (S — 0), with U3 be- 
ing less irrelevant than U2- In this section we perform 
a more elaborate zero-temperature RG analysis of the 
three-body interaction u^. Our aim is to understand how 
a finite imaginary part of the self-energy is generated be- 
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+ 



2 "-<0- 



2 



- 3 




X 

2 -'■<K^- 



FIG. 5: Schematic representations of the flow equations for the impurity propagator (dashed line) and the three-point vertex 
UA4>h (grey circle). Grey squares and black dots denote the four-point vertex ua and the boson three-point coupling gA4>h, 
respectively. Insertions of regulator derivatives and associated combinatorial factors are not shown explicitly. The wiggly line 
is an abbreviation for the difference between the normal and the anomalous boson Green's function [Qp — Fp). Black lines 
without arrows denote either Qp or J-p and we don't show all combinatorial possibilities explicitly. For details see Eqs. (5.9) 
and (5.10). 



yond the kinematic threshold discussed in Sec. Ill A, thus 
we are going to analyze the flow of the impurity propa- 
gator at a fixed external momentum. The two-body in- 
teraction U2 will be set to zero throughout this section. 

We are using a functional RG formulation to derive our 
flow equations.'''^ In order to keep the analysis transpar- 
ent we use simple parametrizations for the propagators 
and vertices, however. Since the coupling between the 
impurity and the bosons is a four point coupling and 
the fRG flow-equations have a one loop structure, we 
employ the standard trick and consider the flow in the 
symmetry broken phase in order to generate a non-trivial 
frequency- and momentum-dependence for the impurity 
propagator while maintaining a simple parametrization 
of the interaction vertices. '^^ This approach is similar to 
the background fleld method used in field-theoretic RG 
approaches. 

Our starting point is to choose a simple truncation of 
the effective action T ^ ip, h* , b] at momentum scale A 
in the symmetry broken phase 



E +C.C.) 



q.k 



J2 r^r^-k^-p^p +T.Ki- ^a',)^« 



5A 

" 4 

q,k,p q 
q,k,p 

>A E "A(f^<;) {b*bq^k1pk + C.C. 
q,k 



(5.1) 



Here is the expectation value of the bosonic field ip, 
the subscript A denotes scale dependent quantities, we 
use the shorthand notation q = {ii}q,q) for sums over 
bosonic Matsubara frequencies and momenta, and we 
dropped the subscript on the three-body coupling = u 
for notational brevity. Note that even though is zero 



at the critical point, the properly rescaled field expecta- 
tion value (f) will be non-zero and is related to the anoma- 
lous dimension of the bosons. 

The effective action in Eq. (5.1) is a straightforward 
generalization of Eq. (2.10) for (5 = and U2 = to 
the symmetry broken phase. The important renormaliza- 
tions appear in the propagators, which we parametrize as 
(we set the velocity of bosonic excitations to unity (c = 1) 
from now on) 



GA^q = -Zt{nl + ci')-rA-gA<l>l (5.2) 

2 

Ba', = zXiilq - «A(fe + ^)+ «7ASgn(r!,) , (5.3) 



where and are wave-function renormalizations and 
sgn denotes the sign function. Note two important mod- 
ifications of the impurity propagator -Ba,?- First, we ex- 
panded the dispersion relation around a fixed momen- 
tum Q = and its corresponding on-shell frequency 
Q^/ {2m) in order to study the flow of the propagator at a 
fixed external momentum. The group velocity v = Q/m 
is allowed to flow, which in turn captures renormaliza- 
tions of the impurity's effective mass. Second, we in- 
clude a damping term 7, which will flow to non-zero val- 
ues above the kinematic threshold. As will become clear 
later, generating a finite imaginary self-energy 7 is the 
only possible way to cure the mass-shell singularity that 
arises from the large phase space for scattering processes 
between the (almost) linearly dispersing impurity mode 
and the linearly dispersing critical bosonic modes above 
the kinematic threshold. 

Also note that we didn't include a chemical potential in 
the impurity propagator. During all our calculations we 
enforce a vanishing density of impurities by neglecting 
the contribution of the impurity propagator's pole (or 
branch cut, if 7 ^ 0) in all diagrams. This procedure is 
tantamount to adjusting a chemical potential such that 
the density of impurities is zero. 

Apart from the impurity propagator also the three- 
body coupling MA in (5.1) will acquire an important non- 
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FIG. 6: Upper panel: RG-flow of the rescaled three-body 
interaction parameters u'^ and «" (shown as soHd blue and 
dashed red line) as well as the non-rescaled imaginary part 73 
(shown as green dash-dotted line) as a function of the flow pa- 
rameter s = logAo/A, obtained by integrating the RG equa- 
tions in appendix B. The initial conditions were chosen such 
that the group-velocity v of the impurity is above the kine- 
matic threshold for creating particle-hole excitations. Note 
that 7j, has been multiplied by a factor 20 for better visi- 
bility. Lower panel: corresponding RG-flow of the rescaled 
impurity velocity Vs- 



analytic frequency dependence during the flow above the 
kinematic threshold, which can be parametrized as 



Here Fj^ [$] is shorthand for the second functional 
derivative of the effective action with respect to the 

fields, i.e. {T^l\^)k<i = 5''T/{6<^k,a.5^q,0) has a 4 x 4 
matrix structure coming from the field index a, P — 
1,...,4 with $fc = {tljl,tljk,bl,bk), besides its usual 
frequency/momentum-space matrix structure. The trace 
in Eq. (5.5) is taken with respect to all momentum, fre- 
quency and field indices. R\ denotes the regulator func- 
tion 



(i?A)fe,, = diag{R%,R%,RU,Rl,)Sk^, , (5.6) 

and we choose an optimized Litim cutoff for both, the 
bosons as well as the impurity'^'' 



A,fc 



ZX(A2-k2)0(A2-k2 



'A,k 



Q 



(A2 - k^)d{A^ - k^) 



(5.7) 
(5.8) 



UAi^q) ^u'j^ + iu'i sgn(r2q) 



(5.4) 



As customary, we neglect derivatives of Z'f^ and wa in the 
scale derivative of the regulator O^Ra. 

The explicit form of our flow equations for all scale- 
dependent parameters is shown in Appendix B. Note that 
the flow equations for the four boson parameters, ta, 5a, 
(f)A and Zf^ are decoupled from the impurity flow equa- 
tions. This is because we only consider the case of a 
vanishing impurity density where all diagrams with in- 
ternal impurity loops evaluate to zero and thus do not 
give a contribution to the flow equations. For this rea- 
son the properties of the bosons at the critical point are 
described by the usual Wilson-Fisher flxed point with 
critical exponents in the universality class of the 3D XY- 
model. 



Flow equations for all scale dependent quantities can 
now be derived in the usual way from the exact flow 
equation for the effective action^^^ 



dA^Am = ^Tr 



1 



Ff [$]+i?A 



9a i? A 



(5.5) 



Without discussing the standard flow equations for 
the bosons much further, we now turn to the impurity 
flow equations, shown in Fig. 5 in a very schematic way. 
The equations for the impurity propagator and the three- 
point vertex are given by 



dABX\ 
9a(ua0a) 



+ -^p) Kp + 2 ^aK E [(^p - -^p)'S,-p Rl^ + {g,-p - T,^p)Bl R\,^ 
p p 

-UA9A<i>A E [2 Ql + 3 Ql^p + 6 QpJ'l + J'l] Rip 

p 

-2 <0A E [(^P - ^P^'^P Kp + (^P - ^pK ^ip 

p 

-SulgAcl^l J2 [2(^P - ^p)'Bp Kp + (^P - •^p)'^P ^A,P 
p 

-2 uUl E [iGp - TpfBl Rl^ + 2{gp - Tp)Bl R^^ 



(5.9) 



(5.10) 



where R — OaR is the scale derivative of the regulator 
and Gp, J'p and Bp denote the regularized normal and 
anomalous boson Green's functions as well as the regu- 
larized impurity Green's function, respectively. 



Qp = 



iGl]p-RtpY-i9AKI'^? 



J' n 



B 



'A,p 



gA0A/2 

^A,p) 



iGi!p-Rt,r-i9AKm' 



R 



A,p 



(5.11) 

(5.12) 
(5.13) 



From the explicit form of the flow equations in App. B 
it is clear that the three-body interaction u is always 
irrelevant. This is mainly due to the screening of the 
three-body interaction by bosonic particle-hole excita- 
tions, giving rise to a contribution to the /3-function of 
the interaction u which is linear in u and negative, inde- 
pendent of the impurity's properties. 

Even though u always flows to zero, a qualitative 
change in the RG-flow occurs if the group velocity v of 
the impurity is above the kinematic threshold for creat- 
ing particle-hole excitations. In perturbation theory this 
happens for w > c = 1, as discussed in Sec. Ill A. In our 
RG procedure the threshold velocity is slightly renormal- 
ized to Uth = {I + + 5*0^/2)^/-^ in units with c = 1, 
where f * , and 0* denote the fixed-point values of the 
rescaled variables (see appendix B). Indeed, for v > Vth 
the RG generates a finite damping term 7 for the impu- 
rity. Typical RG-flows of the impurity parameters in the 
regime v > Vth are shown in Fig. 6. 

We don't find a RG fixed-point that controls the flow 
of 7 above the kinematic threshold. Instead, 7 flows to a 
finite, non-universal value while u flows to zero. In Fig. 7 
we show the value of 7 at the end of the flow (i.e. at A — 
0) as a function of the initial velocity wa=Aoj showing 
the sharp onset of 7 for v > vth- Fromally, generating 
a flow for 7 is intimately linked to the appearance of a 
non-integrable mass-shell singularity above the kinematic 
threshold. 



VI. CONCLUSIONS 

In this paper we have analyzed the problem of a mobile 
impurity atom interacting with bosons at the superfluid- 
Mott insulator critical point in d = 2 spatial dimensions. 
The effective low energy theory describing this systems 
is qualitatively different if three-body interactions are 
present. We calculated the diffusion constant D of the 
impurity as a function of temperature and showed that 
D does not depend on the impurity mass at low tempera- 
tures. Furthermore we showed that the zero-temperature 
properties of the impurity change if its group velocity ex- 
ceeds the kinematic threshold for creating particle-hole 
excitations. This behavior is not controlled by a RG 
fixed-point, however. Instead the imaginary part of the 



0.15- 



y 0.1 




0.05 



FIG. 7: Imaginary part of the self-energy 7s at the end of 
the RG-flow (A — >■ 0) as a function of the impurity's initial 
group velocity t;A=Ao, plotted for two different initial three- 
body interaction strengths uaq = 10 (blue solid line) and 
tJ'Ao ~ 5 (red dashed line). The threshold value of the velocity 
is renormalized from Vth = 1 to Vth = {1 + f* + g,cj)i/2)^^^ in 
units where the boson velocity is set to c = 1 (see text). 



impurity's self-energy flows to a non-universal value while 
the interaction flows to zero. 
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Appendix A: Derivation of the Boltzmann equation 

In this appendix we derive the collision integral for 
the Boltzmann equation in Sec. IIIB as well as the dif- 
fusion equation, which describes the propagation of the 
impurity at long times. Our collision integral will be 
perturbative in U2, but we want to make sure that ef- 
fects to all order in the boson coupling g are accounted 
for, similar to the perturbative calculation in Sec. Ill A. 
This can be done without invoking the Keldysh formal- 
ism by considering an effective model where the impurity 
is coupled to a bath of phonons. As will be shown in the 
following, the spectral density of the phonons needs to be 
chosen such that it matches the spectral function of the 
bosonic current-current correlation function. The colli- 
sion integral for the Boltzmann equation can then be de- 
rived straightforwardly from the effective phonon model 
using Fermi's golden rule. 

We start from the effective action describing the im- 
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purity coupled to a bath of harmonic oscillators 

+ X! "^k.a ^k,ab*q_khq (Al) 

where Wq, denotes the frequency of the phonon mode ^pk.on 
Ck,a are coupling constants and we use the shorthand 
notation k — (Jlk , k) . In order to match this effective 
action to our original action (3.1) we integrate out the 
phonons and obtain 

k,q,q' ^ 

where we have defined the spectral density of the phonons 
via 

J(k, lo)=ttY^-^5{lj-lj^) . (A3) 

On the other hand, if we formally integrate out the 
bosons in our original action (3.1) to second order in U2, 
we get the same result as in Eq. A2 with 

J(k, Lo) = 2ul ImX^o (k, w) . ( A4) 

Now that we have matched the spectral density of the 
effective phonon bath to our original model we can 
straightforwardly employ Fermi's golden rule and obtain 
the collision integral in Eq. (3.8). 

Close to equilibrium we can solve the Boltzmann equa- 
tion (3.6) by first employing a Fourier transform with 
respect to the spatial coordinate x 

9t/k(q) - /k(q) = /[/k(q)] , (A5) 

and making the ansatz 

/k(q) = /^(q) + (k • q)/fe . (A6) 

Here /"(q) ^ nqexp(— /3ek) and fk describes the devi- 
ation of the momentum distribution from equilibrium. 



This is basically a leading order expansion in spherical 
harmonics where the angular dependence of /k(q) on the 
momentum k only enters through the (k • q) term. Note 
that the impurity density distribution is given by 

k k 

Using the ansatz (A6) we can derive a diffusion equation 
from the Boltzmann equation (A5) for arbitrary linear 
collision integrals of the form (3.7) as follows. Taking 
an angular average of the Boltzmann equation (A5) we 
obtain 

2m 27r 

= /[A"(q)] ^ . (A8) 

Moreover, taking an angular average over the Boltzmann 
equation multiplied by k • q we get 

a,Vl?(q) + f^/l?(q) = ^ / 1^ (k-q)/[/k(q)] , (A9) 

where we have used (AS) to replace / with The rhs 
of the above equation can be evaluated straightforwardly 
for a collision integral of the form (3.7) 

/TT in 
^^(k.q)/[/k(q)]= (AlO) 

[^k.k'/fc - ^k',k-^ C0s6lk,k' fk'] ■ 

k' 

Again, using Eq. (AS) to replace / with /° in Eq. (AlO), 
summing over k and neglecting the dff'^ term, we finally 
arrive at the diffusion equation 

dtn^ = -Dq^riq , (All) 

with the diffusion constant specified in Eqs. (3.13) and 
(3.14). 

The explicit expression of the scaling function for the 
diffusion constant in Eq. (3.15) is 



^> '^{xi,X2) = 2tt^xI J J J ^^I^-^oo^y^ {y + 2z - 2y^{y + z)cos0) , y, X2/xi^ x 



(1 + niy)) (l - cos e) e-y + n{y) (l - cos o) 



y + z 

where Kqq is the scaling function of the current correlator 

K^Qiq,uj,A,T)^TKoo{cq/T,Lj/T,A/T) , 



(A12) 



(A13) 
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and n{y) = (e*' — 1) ^. All of the results derived above also apply to the case of three-body interactions if the current 
correlator KQo{q,uj) is replaced by the density correlator n(q, oj) as well as U2 by ^3. 



Appendix B: Explicit form of the RG flow equations 



It has been shown in Sec. V that the flow equations for the boson parameters decouple from the impurity flow 
equations. Indeed, using the rescaled variables 



7 /7V\l/2 <^'A_ ~ ''A 



A(^-i)/ 

and defining the anomalous dimension of the bosons via 

v = - 



9A 



9A 



dlogZf 



rflog A 

the flow equations for the boson parameters are given by (s = log Aq/A) 



47r V32 



15 



dsfs = (2 - 7?) fs + 



4tt V2(l + f,)3/2 



d-^ + V 1 9. 



225 



ds9s = (3 - d - 2?7) 



4^ V32(l + f.)V2+^(-9^^^)j , 
15 ff.^ 



327r(l + f,)5/2 



1 



l+fs 



(Bl) 
(B2) 

(B3) 
(B4) 
(B5) 
(B6) 



In order to bring the flow equations into an analytically tractable form we expanded them in g^(jyi. This basically 
amounts to an expansion in a small anomalous dimension, which works well in d = 2 spatial dimensions. At the 
Wilson-Fisher fixed point the rescaled variables take the values 



^. = ^« 0.035, ~9*~4>l^^^ 



28 

0.049 , n = ~ -0.228 



(B7) 



Note that the value of the anomalous dimension at the critical point 77* = 1/29 agrees pretty well with the expected 
value for a 3D XY-transition, where rj w 0.038.^^ 



For the impurity flow equations we define the rescaled variables and the anomalous dimension of the impurity as 

7A _~ VK . MA 



7A 



VA 

va ~ — r 1 UA ~ I 



Vb 



dlogZi 



and the flow equations take the form 



27r 



Re 



u's9* 



dlog A 



167r 



3 115,g*0j 



(B8) 



(B9) 



(1+^^)5/2 4(l+f^)7/2 



Re 



2tt 



«7s_ 

1/2 {j2Q,3 



3g..^^ 
27r 



Re 



2/ 3i%(a-i7s)^-3us-hw^(4a2-9i7sa-672) 817 



(BIO) 
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IGtt 



3 115 g.^^ 

(l + f,)-V2 4(l+f,)7/2 



Im 



2tt 



lis 



1/2 




~ii',2 I 3i7s(a - «7s)^ - 3-Ds + u^(4a^ - 9i7sa - 67^) 3i7s 



,--,2\3/2 



9,t;, = 



dsjs = (1 - Vb)ls - 



Vs{a- i^sY + ilsvl 1 



< 2(1 + n) 



wja'^ — 17, 
3.9*0^ , 0^ 



_2\3/2 ^Q,3 



2(l+n)5/2 



47r 



Im 



- 173(0; - i7s) 



lis 



(Bll) 
(B12) 
(B13) 



where we have defined — 1 



g^,4>1/'2 and fixed aU boson parameters at their respective critical values. Again 



we have expanded boson loops in g*0j and moreover we have taken Q 00 for the sake of simplicity, which amounts 
to neglecting the small curvature of the impurity dispersion in (5.3). This does not change the results qualitatively. 
We also mention that the term ^ rjijUs cancelled exactly with the ^ contribution in the fiow equation for the 
three-body interaction ma. 

The particular structure with real- and imaginary-parts in the impurity flow equations originates from the non- 
analytic term ~ ijsgnflq in the boson propagator (5.3). After setting the external boson frequencies to zero, all 
diagrams with internal impurity lines depend on the external impurity frequency flq via the damping term ijsgnflq. 
We can extract the contributions of analytic and non-analytic terms by taking limits fig — ?► 0* from above and below 
and adding or subtracting the resulting expressions, which gives rise to the particular form of the impurity flow 
equations above. 
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